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ABSTRACT 

We show that the production cross-sections of 7 and e in p-p coUision, a.^ 



pp^ 



>7 and cTpp^ci 



respectively, may be kinematicaUy equivalent without back to that of the parent pion app^T^- 
So we can obtain straightforwardly app_).e± once we have app^-y with help of the machine data. 
We discuss in detail the relation between them, and present how reliable for (Tpp_j(,± even in the 
TeV region. Based on these studies, we give also the emissivities of 7 and e^ in the Galaxy, 
and present the positron fraction observed at the solar system, the ratio of e+ to [e+ + e~]. It 
is now in a vital open question that it appears to rise as the energy gets higher, > 10 GeV, far 
beyond the expectation with the standard model today. The aim of the present paper is to give 
the production cross-section of e^ firmly even in the TeV region, in order to see quantitatively 
how much deviated from the standard model in the yield of e and its propagation to the Earth. 

Subject headings: cosmic rays — Galaxy: structure — gamma-rays; positrons, cross section 



1. Introduction 

Nowadays the direct observations of antimatters 
such as the antiproton (p) and positron (e^) are 
available with the on board experiments even in 
the higher energy region > 100 GeV, particularly 
by PAMELA (Adriani et al. 2009a, 2009b) and 
AMS-02, while the full result of the latter is not 
yet reported in this stage. 

Since the early reports of the anomaly in the an- 
tiproton flux around 1980 (Golden et et al. 1979; 
Bogomolov et al. 1979; Buffington et al. 1981), 
they have always brought us something new be- 
yond the standard model in both particle-physics 
and the astrophysics, while difficult to induce the 
origin of the anomaly mainly due to the limited ob- 
servations in removing the atmopsheric p's, mod- 
ulation effects, in statistics, and so on. Unfortu- 
nately we do not catch a signal of the anomaly in 
the p spectrum today even with the sophisticated 
modern instruments developed by, for instance 



BESS (Haino et al. 2004) and PAMELA (Adriani 
et al. 2009a) at least up to 100 GeV, which is ex- 
plained by the standard model in the source and 
the subsequent propagation to the Earth within 
the uncertainties in the choice of numerical pa- 
rameters in modeling. 

Ahernatively HEAT (DuVernois et al. 2001) re- 
ported a possibility of the anomaly in the positron 
fiux in 1995 in connection with the WIMPs 
(Weakly Interacting Massive Particles), indicat- 
ing an excess of e^ in the higher energy region as 
compared to the expectation with the standard 
model, and many obsevations for the antimatter 
search have been performed continuously up to 
now with these early studies. 

Recently PAMELA (Adriani et al. 2009a) gives 
the positron fraction, the ratio of positrons (e+'s) 
to electrons (e~'s) plus e+'s, covering the very 
wide energy ranges with 0.1 ^ 150 GeV with high 
statistics. FERMI (Ackermann et al. 2011) also 
gives the positron fraction in 20-150 GeV consis- 



tent with the PAMELA, while not direct observa- 
tion, but the semi-direct one removing the back- 
ground protons in the all positive charged parti- 
cles with the simulation calculations. All these re- 
sults indicate that the positron fraction rises above 
lOGeVupto 100 GeV. 

Under these situations, the AMS-02 is ready to 
make clear the possibility of positron excess with 
higher statistics and unprecedented precision, ex- 
pected to cover the very wide energy ranges, 0.1 
~ 300 GeV. 

With the dvelopments of these observations 
on the antiparticles, the theoretical explana- 
tions for the anomaly in the positron exess have 
been also extensively studied from the various 
points of view, dark matter scenario (Hisano 
et al. 2009a, 2009b; Ishiwata et al. 2009; see 
uittp://pamela. roma2.infn.it/ for a complete list 
of references), pulsars (Biisching et al. 2008; Hop- 
per et al. 2009; Kawanaka et al. 2010), supernova 
explosions happened in the dense gas cloud (Fu- 
jita et al. 2009), reacceleration of secondaries in 
the source environments (Blasi 2009; Mertsch and 
Sarkar 2009; Stawarz et al. 2010), and other mod- 
els (e.g., loka 2010; Kashiyama et al. 2011; Kisaka 
and Kawanaka 2012). 

It has been pointed out, however, that there 
exist many uncertainties in both observations and 
theories for the anomaly, for instance in statistics, 
the energy resolution, the particle identification 
for the former, and artificial assumptions in the 
elementary processes for the latter, such as the 
cross-section of the dark matter annihilation, the 
mechanism for the escape of the pairs from the 
pulsar, and so forth. However, the uncertainties 
in the former would be much reduced soon by the 
AMS-02 as mentioned above, so that the theoret- 
ical interpretations and the parameter sets in the 
assumption must be much more improved without 
worring about the experimental uncertainties. 

The aim of the present paper is not to give the 
explanations for the anomaly indicated by several 
groups today, but to provide definitely the cross- 
sections for the production of 7 and e in p-p col- 
lision and the emissivity of these elements in the 
Galaxy based on the machine data, covering the 
very wide energy ranges in the projectile proton 
even with LHCf (Adriani et al. 2011), 1 GeV - 
20PeV in the laboratory frame (hereafter L.F.). 
Through these studies, we expect definitely how 



much excess appears in the experimental data on 
the positron fraction as compared to the standard 
model for the production and the subsequent prop- 
agation of the cosmic-ray (CR) elements in the 
Galaxy. The theoretical speculations will be re- 
ported elsewhere in the near future with the com- 
ing resuhs of AMS-02. 

One must worry much about the reliability of 
the production cross-section for 7 and e in p- 
p collision for the machine data, among which 
we should remark particularly the energy distri- 
bution, da/dET^, in pion production, which is 
much more important than the angular (pseudo- 
rapidity) distribution, da/drj^^ with 77^=-lntan2, 
as practically we observe the energy spectra of 
photons in the wide energy ranges. 

However, from the experimental point of view 
in the machine data, it is generally very difficult to 
measure the energy distribution of secondary pi- 
ous, particularly in the high energy region within 
the forward cone in the center of mass frame (here- 
after C.M.F.) because of the concentration of sec- 
ondaries in such limited cone, while rather easy 
for the pseudo-rapidity distribution in the central 
region. 

As compared to the charged pious, it is easier 
for the energy estimation of 7's via tt" even in 
the high energy and forward regions, by setting 
the calorimeter with far distance from the vertx 
point, separable enough in the jammed 7's, for 
instance see the LHCf experiments (Sako et al. 
2007; Mase et al. 2012) for the detail. This is 
the reason why we have focussed so much upon 
the data on 7 itself in place of the parent charged 
pious. In fact, LHCf group (Adriani et al. 2012) 
gives us recently unprecedented data on 7's with 
TeV energy in C.M.F. , see Figure 14 in Appendix 
B, and Paper I (Sato et al. 2012) for the detail. 

We have already presented the production-cross 
section of 7's empirically in p-p collision covering 
the energy region Eq = 1 GeV ^ 20 PeV in Paper I, 
and have confirmed that it reproduces nicely var- 
ious kinds of obaservables, pseudo-rapidity, Fey- 
mann variable in addtion to the energy spectrum 
of 7. In the present work, we apply it for the pro- 
duction cross-section of e^ , and point out that the 
production of 7 and e via pion is kinematically 
equivalent, resulting in that we may not need back 
to the parent pion once we construct the produc- 
tion cross-section for 7. 



It would be worth mentioning in advance that 
we separate the secondary mesons into two compo- 
nents phenomenologically, the pionization and the 
isobar ones, the former corresponding to the cen- 
tral collision with the average transverse momen- 
tum of pt « 300 MeV, and the latter to the periph- 
eral one with small momentum transfer, and/or 
with the Pomeron exchange. 

In the present paper, we do not take into ac- 
count the hard process with the large momentum 
transfer mainly effective around the central region 
in the C.M.F. as its cross-section is much smaller 
than that of the soft one with pt « 300 MeV, and 
we are interested only in reproducing the cross- 
section of 7-ray production, particularly in the for- 
ward region. 

The contribution of the isobar component is de- 
termined so that well reproduced are the experi- 
mental data on the energy dependence of the mul- 
tiplicity of mesons (see Section 4.1). In Section 
7, we present the positron fraction based on our 
semi-empirical production cross-section for 7 and 
e , combining it with the standard propagation 
model from the emission site, r, to the solar sys- 
tem (SS), r0. The theoretical implications of the 
result will be discussed elsewhere in the near fu- 
ture with the coming AMS-02 data. 

2. Production cross-section of secondary 
products in pionizaion components 

Before going to the main subject to obtain var- 
ious kinds of cross-section, we touch briefly the fa- 
miliar rapidity-variable 77 corresponding to a par- 
ticle velocity v {= cj3) since it appears very often 
in the present work, well simplifying the procedure 
in derivation of the cross-section we are interested 



m. 



with 



^=^ln^=ln(7+V^^^), (1) 



1 o « 



v/T^^ 



(2) 



For the practical numerical calculations, we had 
better use the second relation in Equation (1), 
since the first expression often gives rise to diver- 
gence for /3 — > 1 in the very high energy region, 
while the first one is much more common in cur- 
rent textbooks than the second. So we define a 



following rapidity- function 77(7), regarding it as a 
function of the Lorentz factor 7 



77(7) = ln(7 + Vt-1), 



(3) 



which appears often in the present work. 

With use of the rapidity ?7s of a secondary part- 
cile s with mass td-s, we have well-known expre- 
sions for the energy and momentum with c ~ 1, 



Es = TOs cosh ?7s ; P^ = tUs sinh rjs , (4) 



with 



d-qs = dEs/Ps 



(5) 



In the following discussion, we call /3 simply as 
a particle velocity in units of c. 

2.1. Decay products coming from pions 

Let us introduce two energy distribution func- 
tions, ^{Eq,E.^) and ip{EQ,€s), each for a par- 
ent pion TT (=7r°, TT ) with energy E'^ [P^^ in mo- 
mentum) in the L.F. and for its decay product s 



(= "7^ 



a ;; 4i,,ib" 



/^ 



with kinetic energy e^ (= Eg — nis) 



respectively, and Eo is kinetic energy of a projec- 
tile proton in the L.F., 



with 



E, 



^Eo,E^)dE^ 
P ' 



l-n1 



2kI 



(6) 



(7) 



where fig is the velocity of s and Kg — tUs /nriT^. 
Namely, two types of decay, 5 = 7, /x, are equiv- 
alent in the decay kinematics but with the mass 
difference only. Note that Equation (7) holds even 
in the case of «„ = with [nig, jig] — > [0, 1] in tt"- 
27 decay. See Equation (13) for the explcit form 
of E^ in such a decay mode. So (y9(_Eo,es) is the 
common function often appearing in the present 
paper. 

2.2. Energy spectrum of 7-rays via tt"- 27 
decay 

In Paper I (Sato et al. 2011), we present ex- 
plicitly the energy spectrum of 7-rays in the L.F., 
ip{Eo, Ej), based on the raw machine data with- 
out referring to tt*^, pointing out that wc do not 



need to worry about the contributions from other 
neutral mesons such as 77, 77', . . . 

Using ip{E(), Ej) given by Equation (6), we have 
a semi-emiprical energy spectrum of 7-rays 



IeI 



N^^{Eo,E^), 



(8) 



explcit form of which is given in Paper I, after the 
integration over the emission angle 9, 






die 



t^c 



with 



■■p Jul- ^0 + (rg 
c^±=/3-i-/3±ie±"o~''- 



c-^o'^dio, (9) 
(10) 



Here, x ~ E^/Eq, uj ~ cos6, and 



770* = ln(2/3,7,Mp/m,o), (11) 

7?; = ln(2/3,7,£;^/m,o), (12) 

and f3c (7c) is the velocity (Lorentz factor) of 
C.M.F. with respect to the L.F., Mp the mass of 
proton, see Appendix A for [Te^Fg] and [0c, C]- 
Note P* = PclcMp in Equation (11) is the mo- 
mentum of projectile (and/or target) proton in the 
C.M.F. . In Equation (10), we have to keep in mind 
W- — —1 for Ey < i77i^oe~(''o +'''=) (77c: rapidity of 
the C.M.F. with rjc = 7;(7c) in Equation (3)). 

Now, while we have used an approximation 
E~ = Ery in Paper I as we have interested in the 
energy region E^ > ttIttO , we use the exact rela- 
tion given by Equation (7) in the present paper. 



EZ 



E^ 



AE^ 



^TT — Eo, 



(13) 



namely we replace E^ with e^ except E^ in uj± 
given by Equation (10) as it comes froms the kine- 
matical constraint in emission angle. 
Another expression of E^ (= e-^) is 



E^ ~ m^a cosh77.y, 
r]^ — \n{2E^/m^o), 



(14) 
(15) 



which is useful for the upper and lower limits of 
pion energy in the electron production as discussed 
in the next section, remembering that a rapidity 
variable 77 in L.F. shifts by 77 — r/c in C.M.F., where 



77c is the rapidity of the latter frame against the 
former one. 

Corresponding to the above replacement, E^ — > 
e-y, we replace x — {E^/Eq) in the integrand of 
Equation (9) by 



Ea 



AxEl'- 



(16) 



so that x^ ~ X for Eq 3> w^o, justifying that the 
results presented in Paper I are applicable enough 
over wide energy region we are interested in but 
in the low energy region Eq < Mp. 

Here we have to take care of the normalization 
constant 0c which was normalized so that the in- 
tegration over X = [0, 1] for Equation (9) leads to 
the average 7-ray multiplicity N^ in Equation (8) 
(see also Equation (10) in Paper I). Instead, we 
renormalize 0r so that 



ip{Eo,E^)dE^^l, 



with 



E^ = - 777^0 e 



Mnl+Va) 



(17) 



(18) 



where r/* is the maximum rapidity of tt" in the 
C.M.F. given by 



Vl = '7(7^); 7^ = /3c cosh77;. 



(19) 



see Equations (3) and (11) for 7/(7*) and 77Q re- 
spectively, and ^* « 77o for Eq > Mp (/3c « 1). 

See Appendix A for the detail of the relation 
between the previous normalization constant 0c 
and the present one after the replacement of x by 

2.3. Energy spectrum of e^ via 7r^-/i^-e^ 
decay 

The electron and/or positron with energy Eg 
are created by the three-body decay from a fully 
polarized muon, /x^ — > e^ + I'di^e) + f'/j(^'/j), the 
production of which is kinematically equivalent to 
7 (see Equation (6)) through the isotropic two- 
body decay as discussed in the last Section 2.2. 

Now, the energy- angular distribution of elec- 
tron (positron) in the muon rest frame is given 

by 



g{q,9) = f{3-2q) 



1 + ^ cos 6* 



nl-2<7 



3-2q 



(20) 



(Commins 1973; Orth and Buffington 1976), 
where q = 2Ee/m^, and ^ = +1 (—1) for elec- 
tron (positron), and 9 is the angle between the 
electron and the spin of muon in the niuon rest 
frame. 

In practice we need the energy distribution of 
e , /(-EjT J -Ee), for a parent pion with fixed energy 
i?7r in the L.F. after the integration over the energy 
and the emission angle of intermediate muon, tak- 
ing the polarization effect into account. The full 
expression of f{Ej^, Eg) is summarized by Dermer 
(1986), while somewhat tedious to evaluate. 

In the present paper, we resummarize it with 
use of a rapidity variable r]^^, introducing a func- 
tion (j), 

f{E^,E,)dE,^^{ri^,qe)^, (21) 

and in Table 1 we present it explicitly with use of 
three parameters, [qe; ?77r, ?7/j], which are given by 

qe = 2Ee/mf,, (22) 

-Ett = JTiTT cosh?77r, (23) 

with rJT, = ri{"/^), and ry^ = 77(7^) in Equation (3). 

Namely, qe is the electron (positron) energy Eg 
in units of half of the muon mass in the L.F., and 
two rapidity variables, rj^ and fy^, correspond to 
that of the parent pion in the L.F. and to that of 
the daughter muon in the pion rest frame respec- 
tively, with 7p = \{m^/m.n + Tn,r/m^) = 1.039, 
leading to 77^ ^ 0.2783 with /3^ =^ 0.2714. See 
also Table 2 for [gi{Q), g2{Q)] with Q = gee^^''", 
gif^M in 0(7y^,gg), and Table 3 for three parame- 
ters [Go; gi,o, 52,0] appearing in [31,32] together 
with their cxplcit numerical values for the electron 
{£_ = +1) and positron (^ = —1). 

In Figure 1(a), we present the numerical value 
of <l){r]Tr,qe) against qeexp[-{r]Tr + fjt^)] for three 
pion kinetic energies, Ett = E'tt — TO,r = 0.01, 0.1, 
and 1 GeV in the case of electron (green curves) 
and positron (red curves). One finds that 0(?7Tr, qe) 
~ (j)(qee~'^^) for e^r > IGeV, namely 4>{rjT^,qe) 
scales in the form of (/)(i?e/£^7r) in the high en- 
ergy region, as naturally expected. 

One should remark that the normalization in 
Equation (21) holds exactly, taking care of the re- 
striction with qe C [g~,g+] in Table 1 

<P{v^,qe)-f^ = l, (24) 



with 



E- = 0; E+ = -TO^e'''+''^ 



(25) 



Once wc have the normalized energy spectrum 
of e coming from the n -n decay in the L.F., 
we obtain straightforwardly the production energy 
spectrum of c in p-p collision (see Equation (6) 
for $) 



dN, 



dEe 
with 



= iV^±/ 'l>{Eo,E^)f{E^,Ee)dE^, (26) 



E^ — niTy cosh rj^ 



(27) 



where ./Vjri (= ./V^ri (Eq)) is the effective multiplic- 
ity of charged pions, depending on the projectile 
proton energy Eq, and 

V^ ^Ve- Vi,; ?7+ = ^; + ije, (28) 

with ?7e = In f/e, and fj* is the maximum rapidity of 
pion in the C.M.F. given by Equation (19), namely 
1]^ is the maximum one in the L.F.. 

Another expression of E~ in Equation (27) is 
given by 



El 



E„ 



4Ee' 



(29) 



which is the same form as Equation (13) for 7-ray 
energy spectrum. 

Now remarking the relation in Equation (7) 
with [E~ ,E:^] = [e^,-Eo] in the case of 7r°-27 de- 
cay, we have 



a7 



-*(So,k; 



1 



p 

and using Equation (21), eventually we obtain 



(30) 



dNei 

dEe 



= A^. 



'p{Eo,e7T)<PHv7T,qe)dr]^, (31) 



where we use a relation dEj^/P-^ = drj.,^, sec Equa- 
tions (19) and (28) for 77^, and e^^ (= E^^ — 111^^) 
is the kinetic energy of pion in the L.F expressed 
with use of vy^r as 



Ett ~ 2TO7r sinh — . 



(32) 



In Figure 1(b), we present the numerical values 
of (j)'' {tIt^ , qe) against geexp[— (yy^r + fj^)] for three 



kinetic energies of the parent pion corresponding 
to Figure 1(a), and again we find tlie scaling na- 
ture in (j)^ (Ee/ETr) similarly as in (f){Ee/ Et^) in the 
high energy region e^r > 1 GcV. 

Note that the energy spectra of two products, 
7 and e , are both expressed with use of the com- 
mon cross-section ip{EQ,e) given by Equation (9) 
with e=e-y,e^, namely the two production cross- 
sections. Equations (8) and (31), are kinemati- 
cally equivalent in the sense that they are linked 
by a model-independent function <j){riT^ , q^) ^ with- 
out referring back to $(i?o,-Eii-), the very model- 
dependent cross-section. 

3. Production cross-section of secondary 
products in isobaric components 

In the low energy region around E^ « 1 — 3 GeV, 
many isobaric components such as Z\(1232), 
Z\(1410), ... as well as the exclusive channel 
pp — > TT+rf (d: deuteron) become effective for sec- 
ondaries in addition to those coming from the 
pionization component mentioned in the last sec- 
tion. 

In this section we consider their contributions 
in an approximate way so that reproduced are the 
experimental data on the production cross-section 
of secondaries, particularly for tt*^ and 7r+, while 
not effective for other mesons, tt^, K , and so on 
(see Figure 3). 

3.1. Energy spectrum of 7-rays via A-7r°- 
27 decay 

We assume that an isobar with mass Ma is pro- 
duced in p-p collision, disintegrating isotropically 
into tt" (-1- p) or tt"*" (-|- n) with the pion energy 
Et^ in the isobar rest frame, and thus we have im- 
mediately 

Et^ = m^TACOshT^A, (33) 



f^A = ln(TAMA/m^), 



with 



TA = Jl-Mp7Mi, 



(34) 
(35) 



note that A'/a must be set with the condition of 
Ett > niT, . One may remark that ta cosh fj\ in 
Equation (33) is the Lorentz factor 7^^ of the pion 
in the isobar frame as presented in Equation (39). 



So, the energy distribution of pions in the L.F. 
is of the uniform one, straightforwardly written as 



<^A{Eo,E^)dE^ = $o(So)dS., 



(36) 



where $0 is determined from the following nor- 
maliztion procedure. 

The pion energy in the L.F. produced through 
the isobar disintegration is bound to the range, 
E~ <E^<E+, with 



E^ = m^ cosh(?7A + ??c ± ?77t 



leading to 



$0(^0) = 



1/(2P.) 



smh{r]*^ + rjc) ' 



(37) 



(38) 



Here 77^ {P^ — m7r/37r77r) is the rapidity (momen- 
tum) of pion in the isobar rest frame, and ?7a the 
rapidity of the isobar in the C.M.F., each given by 

Vtt == v(%)^ % = TA cosh?7A, (39) 

Vk^vhl)'^ 7a = TA cosh 77,*, (40) 



with 



V*s=HV^/rAMA), 



(41) 



see Equation (34) for tja, and y/s = 2Afp7c (center 
of mass energy). Note that rj'^ + rjc in Equation 
(37) is the rapidity of isobar, ryA, in the L.F.. 

Now, similarly as Equation (6), we obtain im- 
mediately the energy distribution of 7-rays com- 
ing from the isobaric components with [E^ , E^] 



leading to 

dN^- 



^a{Eo,E^) 



dE^ 



dE^ 



N^MEo)[Vp - V-f], 



(42) 



(43) 



where N^ is the effective multiplicity of 7-rays 
(via tt") through the A- isobar disintegration, and 
rjp is the rapidity of projectile proton in the L.F., 
and rj-y the minimum rapidity of pion given by 
Equation (15) to produce 7 with the energy of E^. 



3.2. Energy spectrum of e+ via A-7r+-/i+- 
e+ decay 

As the contribution of e^ from tlie isobaric com- 
ponent is negligibly small in comparison with e^ 
around a few GeV, we consider here the energy 
spectrum of e'^ only. 

Similarly to Equation (26), we obtain immedi- 
ately the positron production cross-section coming 



from the isobar decay, A 



■A*" 






Nf^ 



+ 1 <i>A{Eo,E^)f{E^,E,)dE^, 



(44) 

where N^+ is the effective multiplicity of 7r+ 
through the disintegration of A-isobar, and see 
Equation (27) for E^ . 

With use of $o(-Eo) (see Equation (38)) in 
*&a(£'Oj -S'tt), we obtain immediately 






(45) 



see Equations (19) and (28) for ij^. 

Here we have to take care of the integration 
with respect to 77^ in the above equation, which 
must be perfomed separately with the range of 77^ . 
It depends on the positron energy E,, {— 2m ^q^)^ 
see the restriction q^ C [g~,(7+] in Table 1. 

Now, one must note that the energy spectra of 
the secondary products, 7 and e+, coming from 
both the pionization and the isobaric components 
given by Equations (8), (31), (43) and (45), are 
all expressed with use of 0, (l)\ and ^(pdrj^^ alone, 
apart from the basic function ip{EQ,e), where the 
second distribution function 0^ corresponds to the 
differential of the first one with respect to ryjr , and 
the third function to the integral of the first one. 

See next section for the explicit numerical value 
of the parameters appearing above, Ma, N^, and 
so on. 

4. Multiplicity of secondary particles 

In the last section, we present the produc- 
tion cross-sections for 7 and e in p-p collisions, 
[dN^/dEy, dNf,± /dEe] from the pionization com- 
ponents, and [dN^/dE^, dN^/dEe] from the iso- 
baric ones respectively, where the average multi- 
plicity of secondaries, [N^, N^±] and [N^, ^^]j 



appearing in each cross-section respectively, de- 
pend on the energy of projectile proton Eo . These 
numerical values should be determinded from the 
machine data. 

4.1. Parametrization in the empirical for- 
mulae 

In this paper, we slightly modify the formula, 
N^{Eo), on the average multiplicity presented in 
Paper I in order to apply it also for charged com- 
ponents, TT^ and K^, while the contributions of 
charged kaons are not important practically as dis- 
cussed later. 

As mentioned in Paper I, we should regard 
Ns{Eq) as an effective multiplicity, particularly in 
the high energy region, say Eq > TcV, since the low 
energy secondaries are not effective in the galactic 
phenomena. 

We summarize it below 



N,{Eo) = NoA,{Eo)E^ 



0.115 



with 



En — En — _E. 



th' 



(46) 



(47) 




l-e~ 



-(^cA)" 



(48) 

where numerical values of Nq, -E^j^, ec, . . . , are pre- 
sented in Table 4 for individual secondaries, s = 
7, 7r+, . . . , while [61,62] = [4.53GeV, 1.98TeV] ir- 
respective of s. 

For the multiplicity of 7 and 7r+ through the 
disintegration of A-isobar, we assume 



N^iEo)^N^cxp 



IhEq/Ma 
InTTi^ 



2n 



(49) 



where N^ ~ [0.61, 0.15] for s = [7r+,7] respec- 
tively, and E^^ — 0.290GcV (« 2m^), the same 
as the threshold energy of 7r+ production (see Ta- 
ble 4). 

We assume Ma ~ 1.25 GeV, leading to [Ej^, ij,^] 
= [0.279GeV, 2.0] from Equations (33) and (39) re- 
spectively, which is determined so that reproduced 
are the experimental data on the multiplicities of 
7 and 7r+ as presented in Figures 2 and 3, particu- 
larly for 7r+ in the low energy region (green dotted 
curves). 



4.2. Comparison with the experimental 
data on the multiplicity 

Let us compare experimental data with the 
present parameterizations on the production 
cross-sections and the (effective) multiplicity of 
both 7-rays and charged mesons (Stecker 1973; 
Dermer 1986; Antinucci et al. 1973; Sato et al. 
2012). See Paper I for the more complete compar- 
ison between them in addition to the multiplicity, 
where presented were various kinds of physical 
quantities such as pseudo-rapidity, Feynman vari- 
able, 7-ray energy ... in both the C.M.F. and L.F., 
covering very wide energies up to the LHC, Eq = 

1 GeV - 20 PeV. 

In Figure 2, we present simultaneously the pro- 
duction cross-section of 7-rays, cTppiV^, (righthand- 
side axis) and its average multiplicity, N^, 
(lefthand-side axis) against Eq, where Cpp is the 
total inelastic collision cross-section (see Paper I). 
One may note that those coming from the iso- 
baric components are not significant in the case of 
7r° -)■ 27. 

In this figure, we plot additionally two points 
(open and filled red squares symbols) at Eq « 
400 TeV based on the most recent LHCf data with 
yi = 900GeV (Adriani et al. 2012) which is ob- 
tained by fitting our empirical curve in C.M.F. . 
See Equation (18) and Figure 14 in Paper I, and 
also Figure 14 in Appendix B of the present paper 
for the fitting result. 

In Figure 3, we show the production cross- 
section and the average multiplicity for charged 
mesons (Stecker 1973; Dermer 1986; Antinucci et 
al. 1973; Sato et al. 2012), corresponding to Figure 

2 for 7-rays, where we present also the empirical 
curves given by Equations (46) and (49) with use 
of the numerical values in Table 4 in both Figures 
2 and 3. 

Here we present the contributions of 7r+ sepa- 
rately coming from both the pionization and iso- 
baric components, and one finds that the latter is 
effective only around a few GeV with N^ ~ 0.61 
(green dotted curves around a few GcV), while 
not effective for 7r° -^ 27 from A^A-disintegration 
with N^ — 0.15 in Figure 2, so that we do not 
present it explicitly as in Figure 3, but only the 
superposed ones (solid curves). 

From Figure 3, one finds that the contribution 
of K^ are on average of as large as 7% in all 



charged mesons for Eq > 10 GeV and negligible 
below, so that the contribution to e coming from 
K would be of the same magnitude. Thus, in the 
practical numerical calculations, we consider only 
TT^ as the source of secondary e^'^'s, but take 7% 
increase for them due to kaon effect into acount 
finally. 

5. Comparison with other numerical codes 
for the elementary processes 

Nowadays many elaborative codes for the ele- 
mentary processes in the galactic phenomena have 
been extensively developed, and become increas- 
ingly important as the new observational data 
such as FERMI and AMS-02 have become avail- 
able with pretty high quality in both the statis- 
tics and precision. In this section we compare 
our calculations on the production cross-sections 
of 7, e in p-p collision with those expected from 
three calculations by (a) Dermer (1986; analyt- 
ical method), (b) Kamae et al. (2006; simula- 
tional one), and (c) PYTHIA-code (Sjostrand et 
al. 2006). 

5.1. Production cross-section for 7 

In Figures 4(a-c), we present the production 
cross-section on 7's in p-p collision, each with (a) 
Dermer, (b) Kamae, and (c) PYTHIA, for Eq = 
10°, 10\ 10^, 10^, lO'^ GeV, while Dermer presents 
only for the first three, and PYTHIA gives only 
for those except the first energy (not applicable in 
the low energy region < 10 GeV). In these figures 
we plot their results in the form of several discrete 
points alone in order to avoid the complexly for 
the comparison. 

We find that they are more or less compatible 
with each other, while one may claim that there 
exists some deviation with each other. The small 
deviation would be not important in study of, for 
instance, the diffused components, 7 and/or e , 
in the low energy region < 100 GeV, and it must 
be masked by the observational uncertainties in 
practice, even if exists. 

However, it might become serious in the higher 
energy region, say more that 1 TeV. In fact one 
may worry about small (but systematic) deviation 
in the energy spectrum in the higher energy region, 
_Eo = 1, 10 TeV, and it maybe amplified more in 
much higer energy. 



This might be due to the fact that the machine 
data they have based are rather in the low en- 
ergy, < TeV, while the present one covers even 
TeV-7's based on the recent LHCf data, see Figure 
14 in Appendix B as an example, and Paper I for 
the detail in TeV and higher. TeV-7 astronomy 
will be available in the near future, particularly 
after the full observation by CTA-program starts 
around 2020 (CTA Consortium 2012). 

5.2. Production cross-section for e* 

In Figures 5(a-c) and 6(a-c) we present the pro- 
duction cross-sections on positrons and electrons 
respectively corresponding to 7's as presented in 
Figures 4(a-c). Our results seem to be similar to 
those by Dermer in the low energy region but elec- 
trons at -Eq = 1 GeV in Figure 6(a), and we find 
again there exists some dicrepancy between ours 
and others in the high energy region. 

Unfortunately we may have no user-friendly ex- 
perimental data for Eq > 5 x 10^ GeV on the pro- 
duction of secondary electrons (positrons) via ir-fx 
decay, but as emphasized in Section 2 the cross- 
sections for productions of 7 and e^ are equiv- 
alent in kinematics, model-independently linked 
with each other. So we should ask the reliabil- 
ity of the production cross-section of e^ for that 
of 7 which is available enough in current machine 
data. In such a sense, those for e^ must be also 
reliable enough over the wide energy ranges, as we 
have confirmed that our parameterization in the 
production-cross section of 7's is applicable well 
even in the LHC energy region. 

5.3. Spectral-weighted moment in 7-ray 
yield 

In practice, important is the following spectral- 
weighted moment for a given projectile energy Eq 

Z.,{Eo)^ f x^-^cTpp^^iEo,xEo)Eodx, (50) 



where (3 is the index of the differential proton 
energy spectrum with E^ , and note d/dE^ ~ 
Egd/dx. Zj(Eo) must be independent of Eq if 
the Feynman scaling holds exactly. On the other- 
hand, it has been well-known for a long time that 
such a scaling nature is surely broken through the 



study of various cosmic-ray components (electro- 
magnetic and hadronic components) in the atmo- 
sphere (for instance, Konishi et al. 1976). 

Let us present the above moment against the 
projectile energy Eq in the case of /3 = 2.6-3.0 
(index of the proton spectrum in the Galaxy), 
corresponding to 7 (= /? — 1/3) = 2.26-2.66 (in- 
dex of the source proton spectrum), in Figure 
7, where plotted together are other several re- 
sults (QGSJET-II, SIBYLL, and Kamae) sum- 
marized by Kachelriefi and Ostapchencko (2012). 
All results indicate that the scaling nature breaks 
slowly, particularly in the higher energy region > 1 
TeV, but not drastically, which is consistent with 
the past studies for the atmopsheric CR diffusions 
in both electromagnetic and hadronic components. 

6. Emissivities of secondary 7 and e* in 
the GalsLxy 

Once we obtain the production cross-section in 
p-p collision, app^s{Eo,Es) (s = 7,e^), we can 
estimate easily the emissivities of 7 and e^ in the 
Galaxy coming from the interaction between CR- 
proton and the ISM, 



qpp^s{Es,r)^ (51) 

Np{Eo,r)[n{r)vapp^s{Eo,Es)]dEo, (52) 



'E„ 

with 



Tpp^s{Eo,Es) _ dNs , , 

o'pp(-tjo) arjs 



(53) 



where n{r) is the gas density at r, and cjpp is the 
inelastic total collision cross-section (see Paper I). 
The minimum projectile energy -Emin is given 
by the familiar form 



i?it=i?. + ^; (^^7,e±), 



(54) 



which is presented in Equations (13) and (29), giv- 
ing the same form irrespective of s. But in the 
practical calculation, we have to take account of 
the threshold energy effect in the pion production, 

namely Smin = Max[£;th,Emi„] ^'^^^ -^th ~ 2m^. 
In this section we present the emissivities of 7 and 
e at several observation points r in the Galaxy. 



6.1. Normalization of the proton and the 
primary-electron fluxes 

As given by Equation (52), we need two criti- 
cal inputs, the absolute proton density, Np{Eo,r), 
and the ISM gas density, n{r), in addtion to 
the production cross-section, app^g, to estimate 
the eniissivity of the secondary product s at r. 

For the former input, we use here the obser- 
vational data on the absolute proton flux at a 
normalization energy En — 100 GeV, which is free 
from the modulation effect as well as from the ob- 
servational uncertainty in both statisics and the 
experimental precision, while scatters consider- 
ably in the higher energy region > 1 TeV. It would 
be, however, improved greatly very soon after the 
full resuhs by AMS-02. 

Let us present the explcit value of the proton 
density iVOp = iVp(r0, ^„), at {Eo,r) = (En, Vq), 

cN®p = 6.16 X 10"^ (m^^s^^GeV^^), (55) 
based on the current proton flux at _Eo = 100 GeV, 



E, 



2.5; 

° dEr 



= 4.90 X 10^ (m"2g-ig^-iQgV^.5^^ 



(56) 

which is consistent with the past (Derbina et al. 
2005) and the most recent PAMELA data (Adriani 
et al 2011). 

In addition to the absolute proton density, we 
need also the absolute electron density at the nor- 
malization energy £'„ — 100 GeV in the following 



cN'- 



1.82 X 10"-'' (m^^s^^GeV" 



(57) 



based on the current electron flux at Ee = 100 GeV 
mainly focussing on FERMI (Ackermann et al. 
2010), PAMELA (see Figure 9), and the EGG data 
most recently revised by Kobayashi et al. (2012), 

E^--^ ^ 1.45 X 10^ (m-^s-isr-^GeV^), (58) 
dEc 

which is important to see the positron fraction to 
all electrons, c+/(c+ + e~) as seen later. 

Here we have to separate cN^^ into two com- 
ponents, one from (1) the primary electrons, and 
the other from (2) the secondary electrons induced 
by the charged pious via tt"^- fi^- e"^ decays. 



while the contribution of the second components 
is of as large as 10% around 100 GeV. 

With the above normalizations for the proton 
and electron spectra, and if we assume the indices 
of source spectra for two elements are common 
with 7, and given by 



Qp{Ep) — Qo^pEp^, 



Qc{Ec) — Qo,cE, 



-7 

c ' 



(60) 
(61) 



=^n®e = ciV®,+ciV5±, 



(59) 



we can obtain the absolute emissivities of them at 
the source, Qo,p and Qo.o, with help of the solu- 
tion of the transport equations, while depending 
on the source index 7 (see Paper II). We summa- 
rize them in Table 5 for 7 = 2.27, 2.37, and 2.47, 
corresponding to the indices at the SS with /3 (=7 
+ a) = 2.6, 2.7, and 2.8 with a = 1/3 respectively. 

6.2. Emissivities of various kinds of ob- 
servables 

For the ISM gas density in Equation (52), we 
use an empirical one for n{r) presented in Pa- 
per II (Shibata et al. 2011), which is based on 
the GALPROP-code (http://www.gamma.mpe.- 
garching.mpg.de./^aws/aws.html), nowadays the 
most standard one extensively developed by 
Strong and Moskalenko (1998). 

In Figure 8 we present the emissivity of 7 and 
e^ all together at two sites, (r, z) — (a) (r0,O) 
and (b) (r©, Ikpc) with r^ — 8.5 kpc, each corre- 
sponding to that at the SS and at Ikpc in verti- 
cal direction above the galactic plane respectively, 
where we show 7's due to two electromagnetic pro- 
cesses together, the bremsstrahlung (EB-7; sky- 
blue curves) and the inverse compton (IG-7; blue 
curves) induced by the interaction between pri- 
mary electrons and the ISRF(see Paper II). 

From these results, we find the emissitivites of 
the secondary products coming from the proton, 
7 (green curves), e+ (red curves), and e^ (pur- 
ple curves), are all dominant for Ei > 1 GeV 
("i" = 7, e='=) at the SS in comparison with 7's 
coming from the primary electrons due to the EB 
(blue-sky curves) and the IG (blue curves). 

The energies at the maximum emissivities of 
secondary products, 7's and e s, are around 
TO^o/2 w 60-70 MeV for 7, and m^±/3 « 30- 
50 MeV for e^ respectively, as naturally expected 
from the two-body decay in tt*^ — > 27, and the 
three-body decay in /i — > e i^D. 
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7. Positron fraction in all electrons, e + e+ 

Now we give the positron fraction in all elec- 
trons (e"*" + e~), which is obtained by combin- 
ing the emissivity gpp^e±(^o,f) at r presented 
in the last section with the solution of the trans- 
port equation P(rQ, E^;r, Ec), where we need the 
absolute intensities for both protons and the elec- 
trons as mentioned in the last section. See Paper II 
(Shibata et al. 2011) for the procedure to combine 
gpp_,c±(£^c,r-) with P(rQ,E®;r,Ec). 

First in Figure 9, we show the energy spectrum 
for all electrons together with the present calcula- 
tions for three spectrum indices, /3 = 2.6, 2.7, and 
2.8, all normalized at 100 GeV, where we present 
curves with four modulation parametrs, $ = 200, 
400, 600, and 800 MV, using the force field approx- 
imation (Gleeson and Axford 1968). See Derbina 
et al. (2005) for the energy spectrum for hadronic 
components, proton, helium, and so on. 

In Figure 10, we present the current experimen- 
tal data on the positron fraction to all electrons. 



e+/(e^ 



together with our calculations for 



three spectral indices, J3 — 2.&, 2.7, and 2.8, where 
again assumed arc four modulation parameters as 
in Figure 9. 

One may claim that PAMELA data and our 
curves deviate significantly from AMS-01, HEAT, 
etc. at low energies < 5 GeV. The deviation is 
understood in the framework of the charge-sign- 
dependent solar modularion, and its effect is well 
studied with the drift model (Moskalenko et al. 
2002). According to Beischer et al. (2009), they 
find the PAMELA data is consistent enough with 
the pre- PAMELA ones such as AMS-01 (Alcaraz 
et al. 2000; Aguilar et al. 2007), HEAT (Beatty 
et al. 2004), CAPRICE (Boezio et al. 2000), and 
TS93 (Golden et al. 1996), using the GALPLOP 
code with the charge-sign-modulation, while in- 
cluding great uncertinties in the modeling as yet. 

Fortunalely they show explicitly the correlation 
between the positron fraction e+/(e+-|-e~) and the 
positron energy E,, every data points with and 
without the charge-sign-modulation based on the 
GALPROP code. Thus taking the charge-sign- 
dependent modulation into account, we give our 
corrected curves in Figure 11, together with the 
corrected data, where four data, AMS-01, HEAT, 
ACAPRICE, and TS93, are aU combined with 
some statistical weights. Now we find good agree- 



ment each other in the low energy region, while 
we should keep in mind that the correction is 
not model-independent one, including much un- 
certainties in modulation modeling such as helio- 
spheric magnetic field, the tilt of the current sheet, 
and so forth. 

Of course, important is the positron fraction 
in the high energy > 5 GeV as often pointed 
out by past works (for example Moskalenko and 
Strong 1998; Delahaye et al. 2010), and clearly 
the enhancement in the fraction appears, while the 
present result seems to be more significant than 
those expected from GALPROP one in the high 
energy. Uncertainties in both data and calcula- 
tions related to the modulation will be reduced by 
AMS-02. 

8. Discussion and summary 

With use of the refined production cross-section 
for e^ given by Equations (8), (31), (43) and (45), 
we give its emissivity q-p-p^e±{Ef.,r) at r in the 
Galaxy, and the positron fraction, e"'"/(e"'" -I- e~), 
at the SS combining qpp^c± with the propagation 
function, P{rQ,E®;r,Ec)- As reported by the 
previous works, we find again the significant dis- 
crepany between the numerical calculations based 
on the standard model and the observational re- 
sults. 

Now, we must touch upon the contributions 
from the heavy components. He, . . ., Fe as we have 
focused only on the proton. Unfortunately nowa- 
days we have not yet enough data on the nuclei 
effects in the multiple meson productions, while 
various theoretical models are proposed, reproduc- 
ing in thier own way the machine data. 

About this critical uncertainty we have intro- 
duced the so called "enhancement factor" , defined 

by 

esir,Es) ^ qi,ii^s{r,Es)/qpp^s{r,Es), (62) 

with s = 7, e^, p, . . ., where gaii^s is the to- 
tal emissivity of secondary s including all CR 
elements (projectiles) as well as the helium gas 
(targets) in the ISM, while Qpp^s from only p- 
p interactions. As presented in Shibata et al. 
(2000, 2007), the magnitude of Eg is approxi- 
mately around 1.5-1.6, depending very weakly on 
{s;r,Es). Several authors (Stephens & Badhwar 
1981; Gaisser & Schafer 1992; Cavallo & Gould 
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1971; Dermer 1986) also give the enhancement 
factor with nearly the same magnitude, but us- 
ing slightly different definition without the energy 
dependence. This tells us that the choice of the 
model in the nucleus-nucleus interaction does not 
give rise to any serious effect. In practice we as- 
sume £s = 1.53 irrespective of s = 7,e . 

In order to estimate the positron fraction, we 
normalize both proton and primary electron fiuxes 
at 100 GeV, reliable enough even today as the 
statistics is rather rich with small scatterings 
around there as seen from Figure 9. However, 
we should not forget that the spectral shape, i.e., 
the spectral index j3 in both components is quite 
essential in both lower and the higher energy re- 
gions, particularly in the latter, while the change 
of the index in the former may be masked by the 
solar modulation. 

In fact, if we assume another spectral shape (or 
energy-dependent one) in the primary components 
either with softer or harder ones, one must get 
another conclusion on the positron fraction. For 
instance, recent data on the proton and the he- 
lium spectra by ATIC (Wefel et al. 2008), CREAM 
(Ahn et al. 2010) and PAMELA (Adriani et al. 
2011) report all the drastic hard spectrum in the 
high energy region > 100-200 GeV, which must re- 
sult in the hard spectrum of positrons in the higher 
energy region. 

In contrast to these critical indications, we as- 
sume here rather conservative ones with /? = 2.6- 
2.8 independent of the primary energy and the 
species. These uncertainties in the observations 
must be removed very soon by the coming AMS-02 
data as repeatedly reffered, which will bring us not 
only the e^ spectra, but also the unprecedented 
spectra for hadronic componets, proton, helium, 
. . ., as weU as the B/C, lOBe/^Be, etc.. 

In the near future, we will study the origin of 
the anomaly in positron fraction from the thcorit- 
ical point of view, based on the present results as 
well as the unprecedented data by AMS-02, where 
the uncertainties in both observation and the mod- 
eling must be much reduced. 
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APPENDIX A 

Summary of the variables defined in Paper I and 

the numerical value of the renormalization 

constant. 

In this appendix, we summarize several vari- 
ables introduced in Paper I, and present the nu- 
merical values of the renormalized constant Oc^ 
comparing with the previous one. 

First we give below parameters in the L.F. ap- 
pearing in Equation (9), 



re-2(72-l)(Mp/po)sin0, 



r„ 



2{^l-l){l- 13c cost 



(63) 

(64) 



In Paper I, putting x* = E*/Tc = 'i.E*/^ (^s: 
center of mass energy) with Tc — Tc/po^ we assume 
< X* < 1, leading to a following normalization 
constant 6*0 (defined as Oc in Paper I) so that we 
have an average 7-ray multiplicity after the inte- 
gration over the all phase space in the distribution 
function. 



1 _ r ^0(T-e*) 



Ou Jo l + Cr^* 
with C = 0.02, and 



d{cos9*); Tg»=TcSm9*, (65) 



Fi{r)= {x*f{l-x*fc-^-'dx\ (^ = 0,1). 
Jo 

(66) 

Po in Te corresponds to the average transverse 
momentum of 7-rays, pt, given by 



Po Jo l + C^-e 



d{cos0* 



/o 1 + Ct-; 



d{cos9* 



(67) 

and the average transverse momentum pt is de- 
termined from the experimental data, depending 
(very weakly) on Eg {=Eo + Eth), given by 



Two of authors (T. S. and R. Y) deeply ap- 
preciate Research Institute, Aoyama-Gakuin Uni- 
versity for helping and progressing our research 



~ 0.0286 

Pi(Eo) = m^Eo 



1 -cxp -1.156 v£;o 



12 



As presented in Equation (9) in the text, we do 
not use an approximation with x = E^/Eo, but 
the exact one with x-y = x + rn^o / AxEq in the 
L.F. as given by Equation (16). Corresponding to 
the replacement of a; by x^, 0q is replaced by 



1 
with 






dx 



- {l-x^Fef _^ 



fe + Cre 



x± = (m^o/^o)e 



Mvl+vc) 



duj, 
(69) 

(70) 



where uj± is given by Equation (10) and see (19) 
for 77* in x±. 

In Figure 12, we present the previous normaliz- 
tion 6*0 (dotted curve) and the renormalized one 
0c (broken dotted one), where we present the ra- 
tio 0c/ 0Q against Eq (solid curve) together to see 
the effect of the replacement of a; by x-y . One finds 
some difference between them in the energy region 
Ef) < 10 GeV, while not effective in the higher en- 
ergy region with the constant ratio (^ 20% higher 
than the previous one in absolute value). 



APPENDIX B 

Fitting of the production cross-section with the 
machine data. 

We have already compared our empirical cross- 
section in detail in Paper I with the machine data 
in the wide energy ranges, Eq = 1 GeV ~ 20 PeV, 
and find the present parameterization reproduces 
nicely the data even for the energy spectrum by 
LHC. However, as presented in Section 2.2, we use 
the exact minimum energy E~ in 7r° decay given 
by Equation (13), not E^ but e-y, which deforms 
slightly the cross-section in the low energy region 
Ef) < 1 GeV, while negligible in the higher energy 
region. 

So in Figure 13, we give again the 7-ray en- 
ergy spectrum with the revised normalization as 
presented in Appendix A in the case of Eq = 
0.97 GeV, where we present both empirical ones, 
the previous one (dotted curve) and the present 
one (solid one). Thus we find the previous one 
does not reproduce the data for E^ > 0.7 GeV as 
naturally expected, while the present one repro- 
duces well the drop due to the constraint in the 
phase space. 



After Paper I, a new data of LHCf with ^/s = 
900 GeV is reported (Adriani et al. 2012), so that 
we additionally show the fitting result in Figure 
14, reconfirming the present parameterization re- 
produces again excellently the production cross- 
section in the extremely high energy region, in con- 
trast to the fitting in the very low energy region 
in Figure 13. 

In Figure 14 two production cross-sections are 
demonstrated with two emission angles, one with 
^ = 39 /irad (solid square) and the other with 
6 ~ 234 /irad (open square) . One must note 
that our parametrization reproduces surprisingly 
well the data with the set [N^,po] — [30.2-35.2, 
188MeV/c], which are plotted onto Figure 2 in 
the text. 

Anyway our parameterization reporduces the 
experimental data from GeV to PeV with a simple 
form given by Equation (6). 
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Table 1 
Summary of the normalized distribution function of c^-spectrum, (j){r]T,,qe) and its 

DERIVATIVE cj)^ {riT^ , Qe) , RESULTING FROM A CHARGED PION DECAY, TT^ — > ^^ — > e^, WHERE THE MUON IS 

CREATED FULLY POLARIZED, LEFT-HANDED FOR 7r+ DECAY AND RIGHT-HANDED FOR TT^ DECAY. IN THIS 

TABLE, WE USE A PARAMETER, q,. = 2Ee/m^, AND TWO RAPIDITIES, 77^ (=0.278) AND 77^, EACH 

CORRESPONDINT TO THAT OF MUON IN THE PION REST FRAME AND THAT OF PION IN THE L.F., 

RESPECTIVELY. SeE TABLE 2 FOR [gi{Q), g2{Q)] AND [g\{Q) , gliQ)] WITH Q = qee^'^^, e^''*'. 



<P{V^,<le 



<f}HV7v,qe 



Qe C [g^ 



■n^ $'7m 



giiqec^"') ~ gi{<ieO ''") 

52 (gee'''') -52 (gee"'''') 

52(e^'') -32(e-^^) -h5i(e"''^) -5i('7eC-'''') 

52 (e^-)- 32 (gee-"") 





5l(gee'"') + (7l(gee- 


-Vt 


) 


[0, e-('"'+"^)] 






gliqee"") + gliqec' 


-rn 


) 


[^-(v^+v^)^ ^-\v^-fi^\] 






gl{qee'^')+gl{qeO' 


-Vr 


) 


[g-|l,r-^Ml^ g+l'77r-^Kl] 


»7,r < 7)m 


ffKgee-"") 






[Q-\v^-fi^>\^ g+l»7x-^Ml] 


Vtt >7?m 


gliq.e-^') 






^Q+\v^-fl,^\^ g+(»7x+<7M)] 













\f,+ (v.+V^)_ „1 





Table 2 
Summary of [gi{Q), gl{Q)] and [g2(Q), gliQ)] appearing in Table 1 with 
Q = qeC^^^ , e^''*'. See Table 3 for [Gq; 51,0, 32,0]- 



gi{Q) and/or gl{Q) 



52 (Q) and/or gl{Q) 



G{Q) and/or G\Q) 



3i(Q) = »i.oG(X) 
5l(Q) = 5i.oGt(Q) 



52 (Q) =52,0 
gl{Q) = 52,0 



In 



Q G(Q) 



27m Go 

, Gt(0) 



Go 



G[Q) = 0"- ^Q'c-«"- 



Gt(Q) = 2Q'--Q-'e-«'^^ 
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Table 3 
Explicit forms of [Gq; gi.o, 52,0], and their numerical values for electron (^ — 1) and 

POSITRON (^ — —1). 



[Go; 51.0, 52,0] C = l(e-) C = -l(e+) 



2 1 + (2g - 3)/3^ 2^„ 

3 1-/3, 



1.163 -0.570 



51,0 = — 7m(3 - C/3m) 3.745 4.490 



m„ 



52,0 = i — ^r^(? + 5/3m) 6.780 1.027 

6 m, 72 - 1 



Table 4 
Numerical values of coefficients appearing in the multiplicity given by Equation (27) for 



VARIOUS kinds of SECONDARIES S. 



s 


iVo 


£;th(GeV) 


ec(GeV) 


k 


7 


8.80 


0.28 


0.021 


1.0 


n+ 


4.40 


0.29 


0.021 


1.0 


TT^ 


3.80 


0.76 


0.286 


0.5 


K+ 


0.38 


2.50 


0.021 


1.0 


K- 


0.28 


15.0 


0.286 


1.0 
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positron 


iliiCron : E^-m^ 
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i/^O- 






: 10 \k\ 


- 


If ^ 


sS^:::::r:::niz 





: 100 


■ 


■ ,' 


:::iii::: 






■ 


■ »""" 






. 


(a) q^-2i:,/,n^ 


■ 


■ 1 


■^^■,i''''t'''^^ 




^=0.278 


■ 




- 




':::i:s^^^^^^^ 


\ 


■ 


1 ..'' 


J : . i ; 


N 


\^ 



0.2 



0.4 0.6 

9eexp[-(%+7ii)] 



0.8 




0.4 0.6 



Fig. 1. — Numerical values for (a) (jyij^-^^qe) and 
(b) (j)^ {rjTT , Qe) ■ For e7r(= -En — m-7r) larger than 
1 GeV, both (p and (/>'' scales in the form of Ee/Et^, 
as geexp(-7yn) ex E^jE^,. 



18 



Table 5 
Numerical values of the coefficients appearing in the source energy spectrum given by 

Equations (60) and (61). 



source flux 



/? = 2.6(7 = 2.27) /3 = 2.7(7 = 2.37) /3 = 2.8(7 = 2.47) 



Qo,p (m-3s~iGeV-i) 
Qo,e (m-^s-^GeV-i) 



6.909 X 10"^* 
2.513 X 10"^° 



1.092 X 10"" 
4.236 X 10"^'^ 



1.928 X 10"" 
7.110 X 10~^° 



- lO^(mb) 



■s 10 ^ 




- lO^Cmb) 



(for right axisy. 10^ 



. 10^ 



: compiled by Stecker, Dermer 

: previous work without LHCf 

: present work with LHCf 

: previous work from charged hadrons 



10" lOi 10^ 10^ Iff" 105 106 10? 
proton kinetic energy; £o (GeV) 



: Ifll 



10" 




10-' 10^ 10' 10' 
proton Idnetic energy; Eq (GeV) 



Fig. 2. — Average 7-ray multiplicity, iV^, for 
the left axis, and the production cross-section, 
a-mciNj, for the right axis against the proton ki- 
netic energy Eq. See Paper I and the references 
therein for the data. 



Fig. 3. — Same as Figure 2, but for charged 
mesons, tt* and K^. In the case of charged pi- 
ous, we present separately two contributions from 
the pionization and isobaric components with both 
green dotted curves, and superposed ones with 
heavy solid (7r+, K+) and heavy dotted (tt", K~) 
ones. 
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10-2 10-' 10" 10' 

y -ray energy; £^(GeV) 



Fig. 4. — Comparison of the present empir- 
ical production cross-section of 7-rays (heavy 
sohd curves) with three numerical results, (a) 
Dermer (1986), (b) Kamae et al. (2006), and 
PYTHIA-code (Sjostrand et al. 2006), for Eq = 
10°, 10\ 10^, 10^, 104 GeV (squared symbols). 
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Fig. 5. — Same as Figure 4, but for the production 
cross-section of positron. 
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Fig. 7. — Comparison of the numerical vaule 
of the spectral- weighted moment in 7-ray yield, 
Zj{Eo), where our curves (heavy solid ones) 
present those in the five cases of spectral indices, 
/3 = 2.6, 2.7, 2.8, 2.9, and 3.0, while simulational 
ones assume /3 = 2.85. 



Fig. 6. — Same as Figure 4, but for the production 
cross-section of electron. 
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Fig. 8. — Emissivities of 7 (green), e"*" (red), and 
e~ (purple) induced by the collision between CR- 
proton and ISM, taking account of the nucleus ef- 
fect (see Section 8 for the explicit corrections), at 
two observtional sites, (a) SS, and (b) 1 kpc above 
the galactic plane (GP) with r© = 8.5 kpc, where 
we calculate them in the case of /3 — 2.7 and 2.8. 
For the comparison, we present 7's induced by pri- 
mary electron interactions with the ISM (EB) and 
ISRF (IC). 
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10" 10' 10^ 10' 

electron (e'+e*) energy; £j(GeV) 
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Fig. 9. — All electron spectrum (e"^ + e~) from 
current data, together wtih the present calcula- 
tions in the cases of /3 = 2.6, 2.7, and 2.8. Numer- 
ical curves are all normalized to Ee — 100 GeV. 
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Fig. 10. — Positron fraction currently available, 
together with our calculations based on the pro- 
duction cross-sections presented in the present 
work, where we assume three spectral indices, (a) 
13 = 2.6, (b) P = 2.7, and (c) /3 = 2.8. See Adriani 
(2009b) and references therein for data. 
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Fig. 11. — Positron fraction after the correction 
for the charge-sign-dependent solar modulation ef- 
fect according to the GALPROP code based on 
the drift model, both for the data and the calcu- 
lations. 



Fig. 13. — Production cross-section of 7's at 
-Eo = 0.97 GeV, where we present empirical ones 
both from the old one (dotted curve) and the re- 
vised one (heavy solid one). 
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Fig. 12. — Comparison of the previous normal- 
ization constant 0o with the present one Go, with 
use of x-y in place of x. 
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Fig. 14. — Comparison with our parameterization 
for the cross-section with those observed by the 
most recent LHCf expriment at y/s= 900 GeV, 
correspondint to Eq « 400 TeV. 
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